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圆柱 结构 中 的 环 癌 SH 波及 洞 体内 
表面 波 不 存在 性 讨论 


邓 良 玉 , 章 小 杉 , 汝 艳 


(西安 理工 大 学 土木 建筑 工程 学 院 工 程 力学 系 ,710048 西安 ) 


摘 要 :目前 有 关 圆 柱 结构 中 的 环 向 SH 波 研究 仅 限 于 实心 柱 体 与 圆柱 党 体 结 构 。 由 于 无 法 直接 获 
得 圆柱 洞 体 中 传播 的 SH 表面 波 的 解析 解 , 对 比 研究 了 实心 柱 体 和 多 种 不 同类 型 的 圆柱 党 体 中 的 环 
向 SH 波 , 从 波 结构 应 变 能 密度 幅 值 分 布 与 党 体 厚 度 的 关系 讨论 洞 体内 SH 表面 波 的 存在 性 。 建 
立柱 坐标 系 下 均匀 弹性 材料 与 功能 梯度 材料 环 向 SH 波 的 波动 方程 ;分 别 求 得 方程 的 贝 塞 尔 函 数 解 
和 有 需 级 数 渐 近 解 ; 进 一 步 计 算得 到 其 频 散 曲线 、 波 结构 与 应 变 能 密度 幅 值 。 结 果 表 明 :震级 数 方法 
可 用 于 计算 圆柱 型 结构 中 的 变 系 数 波动 方 程 , 且 具有 较 高 的 精度 ;圆柱 型 结构 中 环 向 SH 波 的 能 量 
集中 在 外 表面 或 次 外 表面 ,并 随 着 党 体 厚度 的 增加 ,能 量 集中 现象 更 为 明显 ;通过 应 变 能 密度 幅 值 
分 布 规律 推论 出 圆柱 洞 体内 表面 无 法 传播 环 向 SH 表面 波 。 最 后 ,针对 均 质 结构 和 梯度 结构 ,采用 
反 证 法 证 明了 无 法 得 到 满足 洞 体 内 SH 表面 波 豪 减 条 件 的 解析 解 ,从 而 证 明了 该 推论 。 

关键 词 : 环 向 SH 波 ; 贝 塞 尔 函 数 ; 蝴 级 数 渐 近 解 ; 波 结构 ;应 变 能 密度 
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SH waves in cylindrical structures and discussion on the 
nonexistence of circumferential SH surface wave 
in cylindrical cavities 


DENG Liangyu, CAO Xiaoshan, RU Yan 


( Department of Engineering Mechanics , School of Civil Engineering and Architecture, 


Xi'an University of Technology ,710048 Xi'an, China) 


Abstract; Considering that research of circumferential horizontal shear wave (SH wave) has been 
published on the solid cylinders and cylindrical shell structures, aiming to discuss the existence of SH 
surface waves on the cylindrical cavity. Since it is impossible to directly obtain the analytical solution of 
SH surface wave propagating in cylindrical cavity, we compare the tendency of wave structures , strain 
energy density varying along thickness of SH waves on the surface of cylinders and various cylindrical 


shells. The governing equation of SH wave of homogeneous elastic material and functionally graded 
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material in cylindrical coordinate system is established. The Bessel functions solution and the power series 
asymptotic solution of the governing equations of homogeneous elastic materials and functionally graded 
materials are obtained , respectively. Furthermore ,the dispersion curves,the wave structures and the strain 
energy density are calculated. The results show that the power series method can be employed for solving 
wave governing equations with variable coefficients with high accuracy , the energy of circumferential SH 
waves in cylindrical structures is concentrated on the outer surface or subsurface , and the phenomenon of 
energy concentration is more obvious along thickness. It can be deduced from the distribution of strain 
energy density that the circumferential SH surface wave cannot propagate in the cylindrical cavity. Finally, 
for homogeneous material structure and functionally graded material structure, the inverse method is used 
to prove that any analytical solution cannot satisfy the attenuation condition of SH surface wave in the 
cavity. 


Key words: circumferential SH wave; bessel function ; asymptotic solution of power series method; wave 
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structure ; strain energy density 


随 着 无 损 检测 与 超声 器 件 应 用 需求 的 增加 ,多 
种 材料 .不 同 结构 中 的 表面 波 与 体 波 得 到 广泛 关注 。 
涉及 的 材料 类 型 有 弹性 材料 六 DE RUBORE LE E 
压 磁 材 料 ” 、 压 电 半导体 材料 六 等 ;也 涉及 均匀 材 
Ti 功能 梯度 材料 ;结构 包括 平面 圆柱 中 、 
球体 ”椭圆 柱 "" 、 双 曲 柱 结构 后 等。 对 于 圆柱 过 
体 中 的 SH 波 , 研 究 者 相继 分 析 了 壳 体 径 厚 比 与 弹 
性 常数 对 SH 波 频 散 的 影响 王 ”) 。 不 同 参数 对 SH 
衰 逝 波 :六 的 影响 。 随 着 材料 技术 发 展 ,波动 问题 的 
相关 研究 也 由 均匀 材料 拓展 至 非 均 匀 材 料 。 如 : 非 
均匀 空心 圆柱 中 的 瞬 态 波 "1、 非 均匀 压 电压 磁 板 和 
圆柱 中 的 波 … a 

然而 ,关于 柱 体 结构 中 的 环 向 SH 波 研究 仅 限 
于 实心 圆柱 体 与 空心 圆柱 壳 。 这 些 研 究 并 未 涉及 无 
限 厚 洞 体 中 的 SH 表面 波 传播 , 仅 有 部 分 关于 SH 波 
在 圆 形 孔 洞 散射 的 研究 ,是 方程 不 便于 求解 或 者 
是 这 类 波 就 不 存在 呢 ? 文献 [10] 针 对 椭圆 柱 壳 中 
环 向 SH 波 ,明确 指出 能 量 集 中 在 外 凸 表 面 而 非 内 
凹 表 面 。 如 果 对 于 均匀 壳 体 结构 ,能 量 也 只 能 集中 
在 圆柱 壳 体 的 外 表面 ,那么 对 于 非 均匀 壳 体 ,情况 又 
如 何 呢 ? 考虑 到 功能 梯度 材料 (functionally graded 
material , FGM) 是 非 均 质 材料 的 代表 性 材料 ,所 以 可 
以 考虑 模型 中 包含 功能 梯度 材料 ,那么 其 层 合 结构 
中 的 SH 波 传 播 问题 则 具有 一 定 的 代表 性 。 
功能 梯度 材料 是 一 种 组 分 沿 某 一 方向 连续 变化 
的 非 均匀 材料 ,根据 材料 组 分 的 变化 形式 可 分 为 寡 
律 型 、S 律 型 以 及 指 率 型 。 研 究 功 能 梯度 材料 中 
的 波动 问题 时 , 由 于 波动 控制 方程 的 复杂 性 ,很 难 直 
接 求 得 解析 解 。 目 前 求解 功能 梯度 材料 中 的 波动 问 
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题 ,常用 方法 有 解析 方法 和 数值 方法 。 常 用 的 解析 
方法 有 Legendre 级 数 3:77. Wentzel-Kramers- 
Brillouin 15: 7 FEAR PRU MRA ID Ag 
除了 解析 的 方法 ,数值 解法 也 常用 于 求解 功能 梯度 
材料 中 的 波动 问题 。 目 前 常 有 的 数值 方法 有 传递 矩 
阵 法 "有限 元 方法 “等 。 近 年 来 ,关于 功能 梯 
度 材料 及 其 结构 中 的 导 波 传播 得 到 进一步 发 展 。 
YU 等 名 研究 了 初始 应 力 对 FGM 圆柱 中 轴 向 波 特 
性 的 影响 。WU 等 ' 引 则 对 初始 受 压 的 FGM 圆柱 中 
的 轴 对 称 波 进行 了 研究 ,发 现 FGM 空心 圆柱 的 内 外 
压 差 决定 导 波 的 频 散 。SHEN 等 发现 FGM 圆柱 
壳 中 环 向 SH 波 的 截止 频率 近似 于 等 差 数 列 。 
ZHANG 等 5 采用 具有 更 高 精度 的 模型 计算 了 黏 弹 
TE FGM 圆柱 中 的 周 向 导 波 。 

本 研究 采用 和 窜 级 数 方法 ,研究 圆柱 型 结构 中 环 
[i] SH 波 的 传播 问题 ,基于 弹性 力学 基本 方程 ,得 到 
不 同 圆柱 结构 中 均匀 材料 的 贝 塞 尔 函 数 解 和 功能 梯 
度 材料 的 寡 级 数 渐 近 解 。 研 究 了 圆柱 结构 中 环 向 
SH 波 沿 径 向 的 能 量 分 布 , 从 能 量 角度 分 析 弹 性 波 能 
否 在 无 限 厚 洞 体 表面 传播 的 问题 ,以 期 为 超声 无 损 
检测 提供 理论 基础 。 


1 基本 方程 


图 1 为 几 种 典型 圆柱 型 结构 ,其 中 (a) ~ (AT 
别 表 示 均 匀 圆 柱 体 、 功 能 梯度 圆柱 体 均 勺 圆 柱 壳 、 
功能 梯度 圆柱 充 .内 渐变 型 圆柱 壳 与 外 渐变 型 圆柱 
壳 。 图 中 ,圆柱 体 半 径 为 R; 圆 柱 壳 体 的 内 外 径 分 别 
是 RR 和 RR; 内、 外 渐变 型 圆柱 党 的 层 间 处 半径 为 R; 
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功能 梯度 材料 的 材料 参数 沿 径 向 变化 ,是 坐标 7 的 
函数 。 设 SH 波 沿 环 向 传播 , 柱 坐标 系 下 位 移 函 数 
表示 为 

u, 2 0,u, 2 0,w = w(r,0,t) (1) 
IUP, u, vu 、w 分 别 代 表 径 向 \ 周 向 和 轴 疝 的 位 移 
函数 。 由 于 基本 方程 形式 仅 与 材料 有 关 , 下 面 的 推 
Seb ,功能 梯度 材料 中 的 各 类 物理 量 不 用 上 标 ; 而 均 
xx £- H m 


(a) 均匀 圆柱 体 (b) FGM 圆 柱 体 (c) 均匀 圆柱 过 


(d) FGM 圆柱 过 


(e) 内 渐变 型 圆柱 壳 
图 1 几 种 典型 圆柱 型 结构 
Fig. 1 Several typical cylindrical structures 
柱 坐 标 系 下 ,应 变 分 量 与 位 移 函 数 的 关系 可 由 
几何 方程 得 到 


E 外 渐变 型 圆柱 壳 


Yre 5a Ye (2) 
式 中 , y,, 和 Yo 分 别 表示 r-z 和 0-z 方向 上 的 剪 切 应 


AER ,其余 方向 上 的 剪 切 应 变 和 正 应 变 均 为 堆 
材料 的 本 构 方程 为 


Ta = Mya 5 Tg = MY& (3) 
式 中 , To TS 代表 两 个 剪 切 应 力 分 量 ， 其 他 应 力 分 
量 为 零 ;凡是 前 切 模 量 , 对 于 功能 梯度 材料 , 剪 切 模 


t u 为 半径 > 的 函数 , 即 w = ur 。 
SH 波 在 圆柱 型 结构 中 传播 需 满足 运动 方程 , 即 


ôo, 1 Ty OT, 0O,- du, 
T F —Pp^l2 
ðr r 00 ðz r ðt 
0T, 100, ITa 27, du, 
+ + + 2 4 
ðr r 00 ðz r k ot? (4) 
0T, 1 ÔT o. 0g, Tu Fw 
+ + + =p 
or r 00 0z r ðt 


式 中 ,p 是 材料 密度 ,也 是 半径 7 的 函数 , 即 p =p(7); 
基于 上 述 位 移 假设 ,前 两 个 运动 方程 自动 满足 。 
将 式 (2) 和 式 (3) 代 入 式 (4) 得 到 由 位 移 函 数 表 
示 的 运动 方程 , 即 
du ðw w 


+ 
dr ðr or? 


2 2 
L 9w pow - à w (5) 


r 909? r ór p ot 


报 LIE: 
如 果 材 料 为 均匀 材料 , 则 运动 方程 退化 为 
p 239 de m £n x - n (6) 
研究 圆柱 型 结构 中 的 环 向 SH 波 , 位 移 函 数 除 
需 满 足 式 (5 ) 或 式 (6) 中 所 示 的 运动 方程 ,还 需 满 足 
表面 应 力 自由 边界 条 件 、 层 间 位 移 、 应 力 连 续 边界 条 
件 或 解 的 唯一 性 条 件 。 对 于 不 同 圆柱 型 结构 所 需 满 
足 的 条 件 不 尽 相 同 ,考虑 式 (2) RUSO ,针对 图 1 
中 的 不 同 结构 ,可 分 别 写 出 由 位 移 表 示 的 边界 条 件 、 
连续 性 条 件 或 解 的 唯一 性 条 件 。 
圆柱 体 (图 1a ~b) ,均匀 材料 加 上 标 H 


边界 条 件 TU) ”=0, 叭 一 性 条 件 w |,-。=0 
圆柱 壳 ( 图 1d ~ e) ,均匀 材料 加 上 标 H 


内 渐变 d um 外 渐变 型 圆柱 壳 类 似 


边界 条件 (ae) mou W) uc 
连续 性 条 件 w| = a las | 
(2| | 


2 问题 求解 


对 于 功能 梯度 材料 和 均匀 材料 ,SH 波 环 向 传播 
的 位 移 函 数 试探 解 可 设 为 
WCr,0,D = expliCER9 - oD IWC 
pm = exp[i(kRO — oD1W' co) 
AP: i = V-1; JEDE o 为 频率 ; Woo 和 
WECO 是 待 求 的 位 移 函 数 的 幅 值 。 
将 试探 解 式 (7) 代 入 式 (5) 和 式 (6), 原 偏 微 分 
e im A 


(7) 


人 
(8) 
WARME. 
人 
式 中 :下 是 位 移 ; a = ER; B = pue. 
2.1 控制 方程 中 位 移 函 数 的 解析 解 
2.1.1 均匀 材料 中 位 移 函 数 的 贝 塞 尔 函 数 解 


式 (9) 是 贝 塞 尔 方程 ,其 解 为 


投稿 网 站 :http://cjam. xjtu. edu. cn 微 信 公众 号 :应 用 力学 学 报 


"d 


& 2 期 


W'(r) = CJJ,(gr) +CY (Br) (10) 
式 中 :J ERN SER PR; Y Ti S8 28 Ul ER PR 
数 ; C, 和 C, 是 待定 常数 。 


2.1.2. 功能 梯度 材料 圆柱 体 中 位 移 函 数 的 暴 级 
令 r = r/R , 式 (8) 转 化 为 
2 
(t eu iu * pu TE + po’ - uK)R'W = 0 
dr dr dr 


(11) 
设 材料 参数 可 写 为 震级 数 的 形式 , 即 
u= AIS ,0 = AGI (12) 
n=0 n=0 
如 果 材 料 参 数 变化 规律 确定 ,2 RT b, 0 可 以 
确定 。 
设 式 (11) 可 以 用 寡 级 数 表示 , 即 


W= Yar (13) 


式 中 , m 和 a, 是 待 求解 量 。 
将 式 (12) 和 (13) 代 入 式 (11) 得 到 递 推 方 程 


NC £n - j) (m +n) - (kR)! ]b/ a, P = 


-Ro 2 pg uu (14) 
t 
Sn = 0 得 到 待 求解 量 m 的 方程 , 即 
bi ao En? - CkR?120 (15) 


解 之 得 , m -kR, 
对 比 递 推 方程 中 x 的 系数 ,得 到 决定 a, 的 递 
推 级 数 方程 , 即 


n 


S [(m+n-j)(m +n) - (KR)! ]b/"a, ; + 


j-0 
Ro’ > Das = 0 (16) 

j 
Xia, las di 0;a 为 待定 常数 。 由 于 式 (16) 


P, a, 的 递 推 关系 为 线性 关系 ,因此 可 以 先 对 ao WR 
值 为 1, 然 后 由 递 推 关系 计算 其 他 系数 。 
R3) 可 等 价 为 


W- C, Y a,r“ n+kR> (17) 
n=0 
式 中 , c, 为 待定 常数 。 对 于 其 他 n 值 , a, 由 式 (16) 
确定 。 
2.1.3 功能 梯度 材料 圆柱 党 中 位 移 函 数 罕 级 数 解 


令 h =R-R,r= (r-R)/R-R)-= 
Cr -= ROZh , 式 (8) 转 化 为 
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. Ry dudW /. Ry dW 
— 一 一 十 m RN 
[peu d7 d7 [+ 到 

u kR)? W + (7 + d =- (Fh + Rpa W 


(18) 
设 材 料 参数 可 写 为 震级 数 的 形式 , 即 
u= FOP, p E XXL (19) 
如 果 材 料 参数 变化 规律 确定 , b A 5b,” 可 以 
确定 。 
设 式 (18) 的 解 也 可 用 需 级 数 表示 , 即 
W = Y a,r" (20) 


将 式 (19) 和 式 (20) FICA St C18) rp Sa e t 
数 方程 ,对 比 递 推 级 数 方程 中 的 系数 ,得 到 决定 
a 的 递 推 级 数 方程 


R 
{esr 2) - GRY] tanga Qj 1) 


R 2 
(m € 1) tao (n cj 2) (n D] 


Y PU +w Yo? (ka,j 2hRa, ji + 
j=0 j=0 


Ranj) =0 (21) 
式 中 : Q 440.» 都 为 0; aoa; 为 待定 常数 。 
由 于 式 (21) 中, a, 的 递 推 关系 为 线性 关系 ,为 
了 求解 递 推 级 数 方程 中 的 待定 常数 ,建立 如 下 和 矩阵。 
[ayay] = 了 (22) 
式 中 :j=1,2 E32 x2 的 单位 矩阵 , 则 式 (20) 可 等 
价 为 


W=C $ anr" +C, Ñ apr (23) 
n=0 =0 


xm, C, , C; 是 待定 常数 。 当 n = 0 和 7 = 1 时 ,av 
由 式 (22 ) 确 定 ;对 于 其 他 的 n 值 ,a 由 式 (21) 确 定 。 


2.2 不 同 结构 中 环 向 SH 波 的 频 散 方程 


将 不 同 结构 中 位 移 函 数 的 解析 解 代 入 相应 的 边 
界 条 件 中 ,由 C; 的 系数 矩阵 行列 式 为 0, 即 可 得 出 相 
应 的 频 散 方 程 , 即 波 数 左 与 频率 w 的 关系 。 

2.2.1 均匀 圆柱 体 

均匀 圆柱 体 中 波动 控制 方程 的 解 为 式 (10)。 
第 二 类 贝 塞 尔 函 数 在 r = 0 时 是 发 散 的 ,由 解 的 唯一 
性 条 件 得 到 C，= 0 ,由 此 得 到 


W'(r) = CJ gr) (24) 
RP, a = kR,B = V p'/u"w 。 
得 到 频 散 方程 为 
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Ja (BR) -Jan (BR) =0 (25) Wi = C3J,(B,) + C,Y.(B.) (31) 
对 于 洞 体 情况 ,其 中 一 个 边界 条 件 是 衰减 条 件 ， , T 
式 中 ， em kR, P z /, [9] 
BI reo ,位移 为 零 。 但 由 该 边界 条 件 无 法 得 到 两 。 MuR DOPO Y 
个 待定 常数 M C 的 关系 ,因此 无 法 找到 满足 无 ol ed ed (32) 
穷 远 处 衰减 条 件 的 贝 塞 尔 函 数 解 。 d A : Li í 


2.2.2 功能 梯度 圆柱 体 
功能 梯度 圆柱 体 中 波动 控制 方程 的 解 为 式 
(17) ,由 此 得 到 频 散 方程 , 即 
| (26) 
2.2.3 均匀 圆柱 党 
均匀 圆柱 壳 中 波动 控制 方程 的 解 为 式 (10) , 令 
TIRAM R = R, ,内 径 为 R, ,由 此 得 到 
W'(r) = CJ (Br) + CY (Br) (27) 
RP: a = kR,,B = Vp" "uo o 


则 频 散 方程 为 
ATE (28) 


式 中 
Qu = Joi(BR,) = Jan BR,) 
Qi = Y, 4 C8R,) = Yari (BR.,) 
Qa = Jaa (BR) - J,4 (GR) 
Q5 = Yai (BR) - Yani (BR, ) 
2.2.4 功能 梯度 圆柱 沈 

功能 梯度 圆柱 完 中 波动 控制 方程 的 解 为 式 
(23) ,S EIME R = R, ,内 径 为 R, ,由 此 得 到 

W = Yu x Y aur (29) 


得 到 频 散 方程 
ol 0 Slo eh (30) 


W - C, > aar’ +C, > Qr" 
n=0 n=0 
W'(r) = CJ (Br) + C,Y,(Br) 


oo 
式 中 : Qu 724505 = ap; Qa = > da ;0» = 
n=0 
oo 


am o 
n=0 


2.2.5 内 渐变 型 圆柱 沉 

内 渐变 型 圆柱 壳 中 内 层 波 动 控制 方程 的 解 为 式 
(23) ,外 层 波动 控制 方程 的 解 为 式 (10) , 令 壳 体外 
ft R = R,, 内 径 为 R, , 壳 体 的 层 间 处 半径 为 尺 , 由 此 
得 到 


oc oc 
W = C, Ya, r +C, > an r" 
n=0 n=0 


Qi; 724,05 = a5 
Q4 = Ja GR,) = Jasi (BR,) 
Qa = Yai (BR, ) E Y a (BR,) 


Qi = X = È an 
Q33 Y J.C8R) Qu T Y, (BR) 


Qa EX R nA [3 P 
Qo = RR 


Qu z -AÉ (BR) = Jana (BR)! 


Qua 3 -ElY, (gn) 一 Y 4 (BR)] 


其 他 项 数 为 0。 外 渐变 型 圆柱 壳 的 频 散 方程 与 
内 渐变 型 圆柱 这 的 类 似 , 这 里 不 再 歼 述 。 


3 数值 算 例 及 分 析 


3.1 材料 参数 


选用 陶瓷 和 铬 作为 均匀 弹性 材料 。 陶 瓷 的 材料 

参数 为 pi = 3900 kg/m , u; = 118. 11 GPa ; 铬 的 

材料 参数 为 : p) = 7190 kg/m" ,uy = 102.5 GPa l, 

设 功能 梯度 材料 由 陶瓷 和 铬 复合 而 成 ,材料 的 
体积 分 数 沿 径 向 变化 ,参数 变化 满足 

= gir) + gw(r) 

g? (r) = gir) + gy Cr) 

式 中 : vi +v, =1 ,vs +v, = 1, g” (r) Mg” (r) 分 

别 为 FGM 圆柱 体 和 FGM 圆柱 壳 的 材料 参数 ;8 和 

g2 分 别 为 陶瓷 和 铬 的 材料 参数 ;wm Los 表示 体 


积分 数 。 


(33) 


本 研究 的 数值 算 例 中 设 
m. r-R, 
ANE Li RR 


3.2 ”震级 数 方法 与 解析 解 对 比 
对 于 均匀 圆柱 体 和 圆柱 壳 体 中 环 向 传播 的 SH 
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波 , 均 可 以 得 到 波动 控制 方程 的 解析 解 。 将 解析 解 ”0. 998 m, 外 径 为 1 m 的 均匀 陶 资 圆柱 壳 体 前 两 阶 相 
和 和 震 级 数 渐 近 解 的 频 散 曲线 结果 进行 比较 。 图 2 ”速度 频 散曲 线 。 从 图 2 中 可 以 看 到 ,采用 需 级 数 渐 
(a) 给 出 了 半径 为 1 m 的 均匀 陶瓷 圆柱 体 前 两 阶 。 近 方 法 和 解析 解法 求 得 的 相 速 度 曲 线 结果 相同 ,从 
模 态 的 相 速度 频 散 曲线 。 图 2(b) 给 出 了 内 径 为 。 而 验证 了 知 级 数 渐 近 方法 的 正确 性 。 


20 000 12 000 

解析 解 
.. 16 000 — 10 000 。 TARIE 
E E 
X 12000 z 


8 000 


4 000 


4 8 12 16 20 2 000 4 000 6 000 8 000 10 000 
无 量 纲 波 数 kR 无 量 纲 波 数 kR， 
(a) 圆柱 体 (b) 圆柱 壳 


加 2 过 级 数 法 与 解析 法 结果 比较 
Fig.2 Comparison of results between power series method and analytical method 
| 由 图 3 可 知 :功能 梯度 圆柱 体 的 相 速度 低 于 两 
dlc 个 相应 的 均匀 圆柱 体 的 相 速 度 ;功能 梯度 圆柱 壳 的 相 
计算 圆柱 型 结构 中 环 向 SH 波 前 两 阶 相 速 度 频 “速度 则 位 于 两 个 相应 的 均匀 圆柱 壳 的 相 速 度 之 间 ; 内 
散曲 线 ,圆柱 体 与 圆柱 壳 体 的 尺寸 不 变 , 内 渐变 型 圆 “渐变 型 圆柱 壳 的 相 速 度 与 内 层 陶 盗 . 外 层 铬 的 圆柱 壳 
柱 壳 和 外 渐变 型 圆柱 壳 的 内 外 径 分 别 取 为 0. 996 m ”的 相 速 度 逐 渐 趋 于 一 致 ;外 渐变 型 圆柱 过 的 相 速 度 则 
和 1 m, 层 间 半径 为 0.998 m, 要 大 于 内 层 陶瓷 外 层 馈 的 圆柱 过 的 相 速度 。 


12 000 


10 000 
Ta Ta 
E E gooh 
m S 
E =Æ 6000 
4 000 = 
4 8 12 16 20 2000 4000 6000 $8000 10000 
无 量 纲 波 数 kR 无 量 纲 波 数 KR， 
(a) 圆柱 体 (b) 圆柱 过 
7 000 7 000 
—e—FGM ——FGM 
=- -均匀 材料 --- 均 匀 材 料 
7 m 
E E 
Kk 5 
fur Lr 


2000 | 4000 6000 8000 10000 2000 | 4000 6000 8000 10000 
无 量 纲 波 数 kR， 无 量 纲 波 数 kR， 
(e) 内 渐变 型 圆柱 过 (d) 外 渐变 型 圆柱 壳 


图 3 不 同 材料 圆柱 结构 中 环 向 SH 波 频 散 曲线 对 比 


Fig.3 Comparison ofdispersion curves of circumferential SH wave in cylindrical structures of different materials 
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3.4 位移 幅 值 


下 面 讨论 环 向 SH 波 沿 径 向 的 位 移 幅 值 分 布 。 
圆柱 体 半径 为 1 m; 圆柱 壳 内 径 为 0.95 m, 外 径 为 
1 m; 内 、 外 渐变 圆柱 壳 的 内 径 为 0. 996 m, 外 径 为 
1 m, 层 间 半 径 为 0.998 m。 计 算 结 果 如 图 4 所 示 , 其 
中 ws 为 圆柱 型 结构 外 径 的 位 移 幅 值 。 从 图 4(a) 、 
(b) (d) (e) 可 以 看 出 ,圆柱 型 结构 中 环 向 SH 波 
第 一 阶 模 态 的 位 移 幅 值 沿 径 向 从 内 到 外 逐渐 增 大 ， 
并 且 随 着 波 数 的 增 大 ,内 外 两 侧 的 位 移 幅 值 的 差 值 
也 增 大 。 从 图 4(e) 看 出 , 环 向 SH 波 第 二 阶 模 态 的 


1.0 


—kR-10 
—-kR-20 


模 态 一 


0.8 


WW, 
o 
o 


0.80 0.85 0.90 0.95 1.00 0.97 0.98 


rim 


(a) FGM 圆 柱 体 


模 态 一 10 


0.8 


w^ 
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(b) R=0.95 m，R,=1 m 的 圆柱 党 的 第 一 阶 模 态 


模 态 一 
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位 移 幅 值 则 在 外 表面 附近 达到 极 值 , 随 着 波 数 的 增 
大 , 极 值 越发 向 外 表面 靠近 。 也 就 是 位 移 幅 值 更 集 
中 在 圆柱 型 结构 的 外 表面 或 次 外 表面 。 

为 了 讨论 位 移 幅 值 与 结构 尺寸 的 关系 ,计算 内 
径 为 0.95 m, 外 径 为 1.01 m 的 圆柱 完 环 向 SH 波 的 
位 移 幅 值 分 布 ,结果 如 图 4(f)、(g) 所 示 , 并 与 (b)、 
Ce) 进行 比较 ,发 现 随 着 厚度 的 增加 , 环 向 SH 波 的 
位 移 幅 值 集 中 在 外 表面 或 次 外 表面 的 现象 更 加 明 
显 。 综 上 推 岂 SH 波 的 能 量 集中 在 圆柱 结构 的 外 表 
面 或 次 外 表面 。 


-e R=500 
- -kR,-1 000 


模 态 二 


0.99 1.00 0.95 096 0.97 0.98 0.99 1.00 
rm 
(c) R=0.95m, R=1 m 的 圆柱 沉 的 第 二 阶 模 态 


1.0r -+ kR,-505 模 态 一 
-- kR=1 010 w / 


Us ——kR;-505 FGM 


x EE x 
.07 ETI * 06 一 KR=400 — s oal 移交 
32 oL UU x ——kR-zio > 
0.2 
0.5 Us : 
0 
04 
0.996 0.997 0.998 0.999 1.000 0.996 0.997 0.998 0.999 1.000 0.97 0.98 0.99 1.00 1.01 
rim rim 
(d) 内 渐变 型 圆柱 壳 (e) 外 渐变 型 圆柱 壳 (f) R-0.95 m，R=1.01 m 的 圆柱 壳 的 第 一 阶 模 态 
1.0p -*-kR,-505 模 态 二 E 
---kR,-1 010 , 
si 
三 起 
0:95.— 0.96 :0:97. 098 1099 100 10] 
r/m 
(g) R=0.95 m, R,-1.01 m 的 圆柱 壳 的 第 二 阶 模 态 
图 4 不 同 圆柱 结构 环 向 SH 波 位 移 幅 值 比 
Fig.4 Displacement amplitude ratio of circumferential SH wave in different cylindrical structures 
3. 5 应 变 能 密度 幅 值 W 一 3€ TuYn 十 Ta Yo? (34) 
前 面 从 位 移 幅 值 的 角度 讨论 了 圆柱 型 结构 中 环 可 以 求 得 各 圆柱 型 结构 的 应 变 能 密度 幅 值 , 采 


向 SH 波 的 能 量 分 布 情况 ,得 出 了 能 量 集中 在 结构 
外 表面 或 次 外 表面 的 推断 ,下 面 从 应 变 能 密度 幅 值 
的 角度 出 发 ,分 析 圆 柱 结构 的 能 量 分 布 。 依 据 弹 性 
力学 的 应 变 能 密度 公式 


用 与 计算 位 移 幅 值 时 相同 的 结构 尺寸 。 计 算 结 果 如 
图 5 所 示 , 其 中 Wa 为 加 柱 型 结构 外 径 的 应 变 能 密 
度 幅 值 。 由 图 5(a)、(b)、(d)、(e) 可 知 ,与 位 移 幅 
值 相似 ,圆柱 型 结构 环 向 SH 波 第 一 阶 模 态 的 应 变 
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能 密度 幅 值 集中 在 结构 的 外 表面 , 且 随 着 波 数 的 增 
加 ,这 种 集中 现象 更 为 明显 。 由 图 5(c) 看 出 , 环 问 
SH 波 第 二 阶 模 态 应 变 能 密度 幅 值 在 圆柱 这 的 次 外 
表面 取得 极 值 ,并 且 随 着 波 数 的 增 大 , 极 值 越 靠 近 
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对 比 不 同 厚度 的 圆柱 这 ,由 图 5(b)、(c)、(f)、 
(g) 可 以 看 到 , 随 着 圆柱 壳 体 厚度 的 增加 , 环 向 SH 
波 第 一 阶 和 第 二 阶 模 态 应 变 能 密度 幅 值 向 圆柱 壳 体 
的 外 表面 集中 的 现象 越 明 显 。 所 得 结果 与 前 面 推测 


的 环 向 SH 波 的 能 量 分 布 是 相符 的 。 
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图 5 不 同 圆柱 结构 环 向 SH 波 应 变 能 密度 幅 值 比 


Fig.5 Strain energy density amplitude ratio of circumferential SH wave in different cylindrical structures 


4 解 的 不 存在 性 证 明 


柱 坐 标 系 下 , 环 向 传播 的 SH 波 场 控制 方程 为 
式 (9) ,可 采用 反 证 法 证 明 满 足 无 限 厚 洞 体内 表面 
波 传播 的 解 的 不 存在 。 

假设 解 存在 , 则 需 满足 衰减 条 件 

r— otf, W" 一 0 

采用 坐标 变换 的 方法 , 令 x = 1/r , 场 控制 方程 
(9) 变 为 
a ÊW" 3 dW” 

dx? dx 

显然 , 当 r 一 oo 时 , 即 * 一 0 时, W" 必须 是 小 量 ， 
式 (35 ) 的 前 两 项 与 - on W" 均 是 W" 的 高 阶 小 量 ， 


X 


+ (Gg) -ax W" z 0 (35) 


而 仅 有 BW" 是 W" 的 同 阶 小 量 ,方程 (35 ) 满 足 。 必 
须要 求 B68”= 0 ,显然 不 成 立 ,因此 假设 不 成 立 ,方程 
无 满足 衰减 条 件 的 解 。 
同 理 , 采 用 该 坐标 变换 方法 ,对 于 功能 梯度 材料 
无 限 厚 洞 体 , 控 制 方程 为 式 (8) ,可 以 改写 为 
3/ dW dW |, 
ux [s dz tT) , "em 


* (po? -x uk’ P )W «0 


(36) 

同样 对 比 小 量 的 阶 数 , 可 以 得 到 满足 衰减 条 件 
的 解 不 存在 。 

另外 ,对 于 层 合 结构 的 无 限 厚 洞 体 ,无 论 内 表面 

是 何 种 材料 ,假设 有 SH 环 向 表面 波 在 洞 体内 表面 

传播 ,最 外 层 ( 某 一 厚度 处 至 无 穷 远 处 ) 通常 简化 为 

9 人 匀 材 料 或 功能 梯度 材料 。 依 据 上 述 证 明 过 程 同 样 


投稿 网 站 :http://cjam. xjtu. edu. en ” 微 信 公众 号 :应 用 力学 学 报 


应 用 力学 学 报 


可 以 证 明 无 法 找到 满足 衰减 条 件 的 解 。 因 此 可 以 得 
到 无 限 厚 洞 体内 表面 无 法 传播 环 向 SH 表面 波 。 


5 结 论 


=] 


对 于 功能 梯度 实心 圆柱 体 与 壳 体 结构 中 环 向 
SH 小传 播 问题 , 均 可 以 采用 震级 数 方法 求 得 其 渐 近 
解析 解 ,但 由 于 边界 条 件 不 同 ,两 者 的 震级 数 解 形式 
上 有 差异 。 通 过 与 均匀 材料 相同 结构 中 的 环 向 SH 
流传 播 的 贝 塞 尔 函 数 解 进行 对 比 ,验证 了 考级 数 方 
法 的 精确 性 。 进 一 步 计算 几 种 典型 壳 体 结构 的 波 结 
构 与 应 变 能 密度 ,对比 分析 发 现 : 环 向 SH 波 的 能 量 
集中 在 壳 体 结构 的 外 表面 或 次 外 表面 。 当 壳 体 厚度 
不 变 , 波 数 越 大 ,或 者 波 数 不 变 而 壳 体 厚度 增加 ,能 
量 集 中 现象 越 明 显 。 对 于 无 限 厚 圆柱 洞 体 , 可 以 视 
作 厚 度 无 穷 大 的 壳 体 ,而 洞 体 无 穷 远 处 的 能 量 为 零 ， 
与 能 量 集中 的 外 表面 或 次 外 表面 的 现象 不 符 , 所 以 
推断 圆柱 洞 体 表面 无 法 传播 SH 波 。 最 后 ,从 数学 
上 严格 证 明了 无 法 找到 满足 无 穷 远 处 衰减 条 件 的 控 
制 方程 的 解 ,与 本 研究 能 量 分 析 的 结论 一 致 。 
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